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Abstract 

The capability of current and future measurements at low and high energy e + e~ 
colliders to probe for the existence of anomalous, CP conserving, tvW dipole moment- 
type couplings is examined. At present, constraints on the universality of the tau 
charged and neutral current interactions as well as the shape of the r — ► i energy 
spectrum provide the strongest bounds on such anomalous couplings. The presence 
of these dipole moments are shown to influence, e.g., the extraction of a s {m%.) from r 
decays and can lead to apparent violations of CVC expectations. 
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Although the Standard Model(SM) continues to do an excellent job of describing 
almost all experimental data[|IJ, many believe that new physics must exist beyond the SM 
for a number of reasons and that it cannot be too far away. In addition to searches for new 
particle production or the observation of rare or forbidden decays of well known particles 
at colliders, an alternative probe for new physics is that of precision measurements. One 
possibility is that the repetitive nature of the family structure in the SM may provide some 
insight into what form this new physics might take. In particular the detailed properties of 
the third family of SM fermions may be the most sensitive probes of a new mass scale which 
is beyond the direct reach of current accelerators. To this end it is important to examine 
the nature of the b and t quarks as well as the r lepton with as much precision as possible. 

In the t case one has a particularly clean laboratory in the hunt for new physics 
due to the added advantage associated with the fact that strong interaction effects can 
only arise at higher order, if at all, depending upon which r property is being probed. 
For this reason we focus on the r lepton in the following discussion. One possible form 
of new physics associated with the r is a subtle modification of how the r interacts with 
the gauge bosons of the SM, i.e., the Z, 7 and W. These new interactions, if their scale 
is sufficiently large, can be parameterized as a series of higher dimension, gauge invariant- 
though non-renormalizable-operators involving the r, the SM gauge fields and, perhaps, the 
Higgs boson as well. Neglecting the possibility of leptonic CP violation, the new operators 
of lowest dimension that one can construct take the form of anomalous magnetic moment- 
type interactions. To explicitly construct these operators using only the SM particle content 
requires us to explicitly introduce the scalar Higgs field M in order to produce the required 
helicity flip while maintaining gauge invariance. Thus we find that anomalous dipole moment 
operators for the r are necessarily of dimension- six. 

In the case of the Z and photon, our normalization for these operators take the 
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conventional form 
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t 7m (v t - a T 7 5 ) + 



— <Vg"(«? - *«f 75) tZ", 



(1) 



2c. 



where as usual a T = —1/2 and v T = — 1/2 + 2sin 2 6 l , J) . Similarly for the photon we can write 
an almost identical interaction structure: 



where Q T = —1. In both cases k(k) corresponds to an anomalous magnetic (electric) dipole 
complex form factor. Of course as is well-known, radiative corrections in the SM can easily 
induce dipole moment-type interactions. In this specific case, since the real part of the 
R's are intrinsically CP violating, they remain zero to several loops whereas the k's receive 



to the absorptive parts of the loop diagrams. In the present paper we will be interested in 
anomalous magnetic dipole type couplings that are over and above those of the SM and are 
comparable or perhaps somewhat larger in magnitude. We might expect that if very high 
mass scales are inducing these anomalous couplings then the parameters k and k arising 
from this new physics will be real and not to be very scale dependent, i.e., their values at 
q 2 = and M§ will be little different. 

Experiments designed to directly probe the couplings kJ ,z and kj ,z have been per- 
formed at LEP and elsewhere 0] with only negative results. In fact, in the Z boson case we 
can use the Z — > r + r _ decay width, the r forward-backward asymmetry and the angular 
distribution of the r polarization to constrain both k z , which we assume to be real, and 
\k z I through a 'radiative corrections' analysis^. Using the data as presented at Moriond 



1997|1|, a modified version of ZFITTER5.O0 and the input values of a s [M z ) = 0.118 and 
a EM = 128.8960 we obtain the 95% CL regions shown in Fig. 1. From this analysis we 





complex, order a, contributions f3[. In both cases the imaginary parts of k and k arise due 
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can conclude that if the anomalous couplings are the only source of new physics the Z-pole 
data tells us that \nf,kf\ < 0.0023. Note that this value ~ a /it, the typical size of a SM 
loop correction. The direct search results yield comparable limits, particularly in the case 
of the CP-violating couplings. The corresponding bounds on k? t and obtained through 
direct means are somewhat weaker. 

The purpose of this paper to examine the analogous CP conserving charged current 
dipole moment interactions of the r which have not been as extensively discussed in the 
literature ||. To this end, we modify the charged current ruW interaction to now be of the 
form 
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P L v T W» + h.c. , (3) 



which parallels the above coupling structures for the photon and Z. Here Pl is the left- 
handed projection operator and is assumed to be real. If one allows for the possible 
existence of /tf' 7 then must also exist with a comparable magnitude due to the demand 
that the original non-renormalizable operators be gauge invariant. This is easily seen in a 
number of ways but perhaps the simplest is to examine a particular realization for these 
non-renormalizable operators as presented in Ref.0: 

C new = ^la„ v [a rB B^ + a rW T a Wr] t r $ + h.c. , (4) 

Here, L is the left-handed doublet containing the r, A is a large mass scale, a T w,TB are a 
pair of parameters, W£ u and are the SU(2)l and U(l) field strength tensors, T a are the 
SU(2)l generators, and $ is the Higgs field which we can ultimately replace by its vacuum 
expectation value. Re-writing these operators in the mass eigenstate basis and employing 
the normalizations of Eqs. (1-3) we see that we can immediately write a pair of extremely 
simple relations amongst the three anomalous magnetic moment type couplings which do 





Figure 1: 95% CL allowed regions for the ttZ anomalous dipole moment couplings. In (a) 
m H = 300 GeV has been assumed and the three curves correspond to m t = 169 (dotted), 
175(solid) and 181(dashed) GeV, respectively. In (b) m t = 175 GeV is assumed for 
m#=60(dotted), 300(solid) or lOOO(dashed) GeV, respectively. 
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not depend upon the scale A: 




-1 tx 



SC t + X ' 




1 te + 1 

2sc t + x ' 



(5) 



where t = sva.6 w / cos9 w = s/c and x = a T w/ct T B- A plot of these ratios is shown in Fig. 
2 where we see that they are typically of order unity except near the points x = — t and 
x — 0. These arguments possibly suggest that the relevant interesting range of \k^\ may 
not be much larger than about 0.001-0.01 since the corresponding anomalous ref couplings 
were found to to be quite small. Again, this range is not much different than what we might 
expect from a typical loop correction of order a/n in the SM or simple extensions thereof. 
Since we will concentrate solely on below, we will henceforth employ the simple notation 



The first point we need to address is what are the current direct limits on k and 
how far are they away from the suggestive range of 0.001-0.01? The most obvious set of 
constraints arise from considerations of lepton universality since the r now interacts with 
the W in a away that differs from the other leptons, which we treat as having only SM 
interactions. The two best places to test r universality are the r lifetime itself and the 
ratio of W decay widths R = T(W — > tv)/T{W — > nv,ev). It is obvious that ratios 
such as T(r ^ Pv)/Y{P — > fj,v), where P is a pseudoscalar meson, will not yield additional 
constraints. Since the matrix element of the hadronic weak current between the pseudoscalar 
and the vacuum is proportional to the momentum of the virtual W, the a^q" piece of the r's 
leptonic current will drop out, i.e., only the longitudinal part of the virtual W contributes. 
Any deviation from universality observed in such ratios cannot be attributed to an anomalous 
moment for the r. 
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For the leptonic decay of the r we find, neglecting the mass of the lepton in the final 



state 



K K 



t(t^Ivv) = r (i + - + 



2 10 7 

where T is the conventional SM result, while for the W —>■ tv decay width we obtain 



(6) 



T(W —> tv) = 7=^(1 - r) 
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r 3k K 2 (l + 2r) 

H 1 - 

2 4 8r 



(7) 



where r = (m T /Mvi/) 2 . Here we see that in the W decay case the k 2 term is dramatically 
enhanced kinematically. It is easy to include in these equations the effect of a CP violating 
electric dipole moment type term, R, by the replacement of the k 2 in the last term of both 
expressions by the combination k 2 + k 2 . k only appears quadratically since neither of these 
quantities are CP violating observables. Using the most recent experimental data|| we are 
able to obtain the following 95% CL bounds, neglecting potential R contributions: 



BpTr/r^ : -2.24 • 1(T 2 < k < 2.31 • 10^ 
R : -2.54 ■ 1(T 2 < k < 2.83 • 1(T 2 



(8) 



Note that the two constraints yield comparable limits implying that we must have |k| < 
0.0283 ~ 0.03, which is somewhat larger that the bounds we obtained for the Z case. We 
chose not to combine these two results since, in principle, k(0), which is probed by the r 
width may differ from K(M^r), probed in W decays. Note that the W decay expressions can 
also be used to obtain a bound on \R\ of a comparable size, but any limit obtainable from 
the r lifetime will be rather poor. Clearly, we need to improve our sensitivity to nonzero 
values of k by roughly a factor of a few to an order of magnitude. The statistics for doing 
so will be available at future runs of the Tevatron and at future B and r/c factories. 



Where else might we be able to probe values of k of order 0.01 or less? In addition to 



modifying the total r leptonic decay width, the presence of k ^ produces a distortion in 
the final state lepton spectrum which in general cannot be expressed as shifts in the Michel 
parameters. We recall that the general nature of the Michel spectrum assumes the absence 
of derivative couplings-something we now have due to the anomalous moment. We find that 
the invariant double differential decay distribution for r leptonic decays in the absence of 
the polarization dependent terms is given by 



where p(k) is the r's(final state lepton's) four momentum and terms of order m 2 /m 2 have 
been neglected. Note that in the k 2 term cubic powers of p ■ k appear which are clearly 
incompatible with the standard Michel spectrum. As above, we can include the CP violating 
electric dipole contribution by the standard replacement k 2 — > k 2 + k 2 , which implies that 
the lepton energy spectrum is not very sensitive to k non-zero as we might expect. There 
are two places where precision measurements of the r lepton spectrum can be made: at 
high energies, sitting on the Z, or at low energies at a B or r/c factory. On the Z an extra 
advantage is obtained due to the fact that the r is naturally polarized. Beam polarization can 
greatly enhance this added sensitivity as has been exploited by the SLD Collaboration |TT|. 
At the Z, the normalized lepton energy spectrum for the r to leading order in k, which seems 
a reasonable approximation since k is small, can be written in the absence of m 2 T jM\ and 
mj/m 2 corrections as 



where z = Eg/E T , PL* is the production angular-dependent effective polarization of the 



dT 



3m 2 T (p ■ k) - 4(p • k) 2 (l - + l^(p ■ k){m 2 T - 2p ■ kf , (9) 

2 2 mi 



dEid cos 6 T £ 




(10) 



r including the effects of the initial e beam polarization as given in Ref.|10[] and f,g are 
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kinematic functions: 

/ = (3 - 3^ 2 + -z 3 ) + «(-- + -z 2 - -z 3 ) , 

5 = + 3^ 2 - |^ 3 ) + + 2^ - |^ 2 ) , (11) 

The behaviour of both f(z) and g(z) due to variations in k is shown in Fig. 3. There are 
several things to observe: first, g is more sensitive to variations in k than is / so that having 
a handle on PJ^ through modifications of the beam polarization is very useful in obtaining 

a greater sensitivity. Second, the k dependent term in / has the same structure as that due 



to a shift in I • p, where p is one of the conventional Michel parameters ||1 If. The k dependent 



term in g is not of the same form as a shift in the corresponding 8 parameter. 

The current bounds on k from the SLD Collaboration are rather poor due to the low 
statistics available. To get an idea of what future constraints on k might be obtainable, we 
perform a toy Monte Carlo study assuming a sample of 20,000 r — > i decays, corresponding 
to about 1.5 million Z's, assuming a beam polarization of |P| = 77.2% and a cut on the 



angular acceptance of | cos^| < 0.75. We follow the same approach as SLD[10| and consider 
four distinct cases depending on whether P is positive or negative and whether the negatively 
charged r goes into the forward or backward hemisphere of the detector. In each case, we 
divide the z range into 10 equal bins, generating data weighted by the statistical errors 
only. Fig. 4 shows the resulting Monte Carlo generated distributions in comparison with the 
expectations of the SM. 

In order to obtain our k constraint, we perform a % 2 fit to the above distributions, 
which were generated assuming k = 0, by using the ^-dependent functional forms above. 
We obtain a best fit of k = 0.06 ± 0.12 at 95% CL with a x 2 /d.o.f. of 41.8/40; this result is 
explicitly shown in Fig. 5. It is clear from this analysis that Z pole measurements will never 
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Figure 3: The kinematic functions f(z) and g(z) for k values of -0.05(dot), O(solid), or 
0.05(dash). 
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Figure 4: Comparison of SM expectations and Monte Carlo data for the lepton energy 
spectrum in r decays at the Z. The top figure corresponds to the sum of the 2 cases where 
the negative r is forward and P > and where the r is backward and P < 0. The bottom 
plot shows the sum of the other two possibilities. In both figures the SM expectations 
are given by the dash-dotted histogram whereas the points represent the Monte Carlo data 
generated assuming the SM is correct. 
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be able to achieve the level of sensitivity we require to probe \k\ values below 0.03. 




K 



Figure 5: x 2 fit to the Monte Carlo data using the ^-dependent form of the leptonic energy 
spectrum in r decays on the Z pole. 



What happens at a B or r — c factory? Here the r is either at rest or is relatively 
slow and a pseudo-rest frame can be defined. In the r rest frame, neglecting the I mass and 
polarization dependent terms, one finds 



ax L 



(12) 



where x = 2Ee/m T and r was introduced above. (As above, the contribution of k can be 
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included by the replacement k 2 — > k 2 + k 2 .) Neglecting terms of order k 2 again we see that 
the normalized lepton energy spectrum in the r rest frame can be written as 

I^ = x 2 [6-4x + /t(4x-3)] , (13) 
r ax 

As in the Z case, we see that by comparing with the conventional Michel spectrum|J the 
effect of k on the normalized spectrum to this order is the same as 8/35p[0]. The current 



best single measurement of p using the final state lepton spectrum comes from CLEO|p~2 
assuming e—p universality they obtain 5p = p— 0.750 = — 0.015±0.015 or k = — 0.040±0.040. 
This is already far better than what we obtained earlier on the Z peak using our toy Monte 
Carlo. Using the current world average value Sp = p — 0.750 = —0.009 ± 0.014|13[ gives a 
comparable result. It is clear that future .B-factories, with more than an order of magnitude 
increase in statistics should begin to probe k values of order 0.01. Fig. 6 shows a plot of the 
function h{x) = f^E for different values of k. Note that the most significant deviation from 
the SM occurs at large values x where the most statistics are available. 

In order to get an estimate of the future sensitivity of low energy measurements to 
non-zero values of k, we return to our Monte Carlo approach assuming a sample of 10 s r 
pairs. Data on the r's lepton energy spectrum is then generated assuming the SM is correct 
and is put into 10 equal sized x bins as is shown in Fig. 7. Only statistical errors were 
included. 

A /t-dependent fit is the performed with the results shown in Fig. 8. Here we see that 
the fitted value of n is determined to be k = (-0.49to;fi) • 10~ 3 at 95% CL with a x 2 /d.o.f. of 
11.9/10. A parallel analysis using a data sample of the same size but with ten times as many 
bins obtained « = (-0.04 ± 0.79) ■ 10" 3 at 95% CL with a x 2 '/d.o.f. of 90.5/100. From this 
analysis it is clear that very large r data samples of this magnitude will allow the probing 
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Figure 7: Comparison of SM expectations and Monte Carlo data for the lepton energy 
spectrum in r decays at a B or r — c factory. As above, the SM expectations are given 
by the dash-dotted histogram whereas the points represent the Monte Carlo data generated 
assuming the SM is correct. 
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re's in the 0.001-0.01 range even if significant systematic errors are present. 
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17.5 

X 

15.0 
12.5 
10.0 

-0.002 -0.001 0.001 

K 

Figure 8: \ 2 fit to the Monte Carlo data using the re-dependent form of the leptonic energy 
spectrum in r decays at a B or r — c factory. 

Are there additional constraints that we can obtain on re from the semileptonic decay 
modes? As we saw above, r decays to pseudoscalar mesons do not provide any additional 
sensitivity to re ^ since the purely transverse <J^q u terms decouple in this case. We thus 
turn to decays of the type r — > Vu where V is either a vector or axial-vector meson or where 
V represents the hadronic continuum at and above the two pion threshold. As is well known, 
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the ratio of the semileptonic to leptonic widths of the r can be symbolically written as 

r(r -» hadrons) f m l s p{s)ds 
1 (r — > ivv) Jo m; 

where p(s) = Imll(s)/ir is the spectral density defined in terms of the appropriate charged 
current correlator II(g 2 ) and R is the kinematic kernel which now depends on k. Defining 
x = s/m 2 we find that R[x) is given to all orders in k by 



1 + 2x - 3kx + - a k 2 x{2 + x) 
R(x) = (1 - x) 2 - f— i— n . ( 15 ) 



1 + 2 « + iQ K 2 

As usual, we can incorporate k contributions by the now standard replacement k 2 — > k 2 + 
in this expression; again, the numerical contributions from the k terms cannot be large. 
The above expression for R(x), expanded to linear order in k and neglecting possible k 
contributions is shown in Fig. 9 and displays reasonable sensitivity to k. In principle, a non- 
zero value of k can be responsible for an apparent breakdown in the usual CVC predictions 
in that the r's weak coupling to the W is no longer solely a mixture of conventional vector 
and axial-vector currents. Of course, the matrix element of the weak hadronic charged 
current is still simply related to the corresponding electromagnetic current by an isospin 
rotation. However, when it is contracted with the r's charged current matrix element, new 
ft-dependent terms arise which are not present in e + e~ — > hadrons. It is clear from Fig. 9 
that positive (negative) values of k will produce a decrease(increase) in the predicted value 
of the relative hadronic branching fraction of the r, R T . As an estimate of the current 
sensitivity we compare the experimental rate for r — > 1-nv with that expected from CVC 



using the low-energy e + e~ data as input [14 in the region of the p resonance allowing for k 
to be non-zero; we find that k = —0.045 ±0.034 at la. Also, crudely, we find that a value of 
k = 0.03(— 0.03) leads to decrease(increase) in the continuum r hadronic width of roughly 



—3(3)%, a value not too much smaller than the leading inclusive QCD correction to the tree- 
level SM electroweak result of a s (ml)/iT. k ^ could thus have important implications to 



the problems associated with a clear understanding of the running of the QCD coupling [15] 
if the value extracted from hadronic r decays 0] is indeed shifted due to these anomalous 
couplings. 
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Figure 9: The kinematic function R(x) as defined in the text assuming K—-0.05 (dotted), 
O(solid), and 0.05(dashed). 



Lastly, we note that a non-zero value for k also leads to a modification of the transverse 
and longitudinal polarization fractions for the C V in the final state of this semileptonic decay. 
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Unfortunately, the k ^ contribution is quite small particularly in the low invariant mass 
region. With x now denoting the ratio m v /m^ we find that to lowest order in k 

(2 - 3*Qs 
T ' 1 + (2 - 3k)x ' 

f l = (16) 

1 + (2-3k)x' V ; 

where is the transverse (longitudinal) fraction of 'V. For \k\ < 0.03 we see that this 

effect is below the 1% level; there is little sensitivity here to anomalous couplings. 

In this paper we have considered how existing and future measurements can probe 
the CP conserving anomalous magnetic dipole moment of the r, k, in charged current in- 
teractions. Searches for the corresponding coupling structure in electromagnetic and weak 
neutral current interactions of r's have been undertaken for some time and interesting limits 
have been obtained. We have argued that if 7 and/or Z anomalous couplings are present for 
the r at some level then the analogous charged current couplings must also be present with 
a comparable magnitude due to gauge invariance. In almost all cases, except for the decay 
W — > ru, the contributions to the observables that we consider from potential CP violating 
terms arising from the electric dipole moment k, if they are at all present, are shown to be 
small. 

We observed that the current best limits on k are indirect and arise from universality 
tests, specifically, considerations of the r lifetime and the W — > tu branching fraction. 
These two independent measurements tell us that \k\ < 0.028 at 95% CL. Clearly, at both B 
and r — c factories, universality tests involving the r should improve these indirect limits by 
factors of order a few to as much as an order of magnitude depending on the size of systematic 
errors. More interestingly, direct experiment sensitivity to non-zero values of k arise from a 
number of sources, in particular, the shape of final state lepton energy spectrum in r — > tvv. 
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We showed that the shift in the unnormalized spectrum cannot be accommodated through 
variations in the Michel parameters. However, to first order in k, the normalized spectrum 
modification was found to be proportional to 5p. A toy Monte Carlo study showed that the 
sensitivity obtainable from future data on a sample of more than 10 6 polarized Z's at SLD 
cannot probe k values at the requisite 0.01 level. A measurement of the same spectrum at 
B and/or r — c factories where 10 8 r pairs are available was shown by the same toy Monte 
Carlo approach to be sensitive to k at the 0.001 level. 

Hadronic r decays were shown to display somewhat more subtle sensitivity to finite 
k and, in particular, decays to a single pseudoscalar were shown to have no sensitivity 
whatsoever since only the couplings of the virtual longitudinal W are being probed there. 
Decays to 'vector' or 'axial- vector' meson final states(ie., the hadronic continuum above the 
2n threshold) were shown to display three-fold k sensitivity via modifications in the overall 
decay rates as well as in the associated invariant mass distributions and the polarization 
of the final state. In particular we showed that the deviations in the value of R T due to a 
non-zero k can lead to an incorrect extraction of the strong coupling constant, a s (m%.), with 
the obvious implications elsewhere. We also found that the anomalous moment terms can 
lead to potentially large apparent violations of CVC. 
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